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Abstract 

We investigate the angular normal modes for small oscillations about an equilibrium of a single- 
component coulomb cluster confined by a radially symmetric external potential to a circle. The 
dynamical matrix for this system is a Laplacian symmetrically circulant matrix and this result 
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leads to an analytic solution for the eigenfrequencies of the angular normal modes. We also show 
the limiting dependence of the largest eigenfrequency for large numbers of particles. 
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I. INTRODUCTION 



We present in this paper an analytic solution for finding the angular normal mode frequen- 
cies of small oscillations about an equilibrium position for a system of interacting charged 
particles confined by a radially symmetric external potential. Experimental applications of 
this analysis might include the following types of systems: laser cooled ions in a Penning 
ions in a cylindrically symmetric mass spectrometer^^, clusters of electrons on the 
surface of liquid heliurrP, and atoms confined in a chemical ring^. We limit our investi- 
gation to two-dimensions, similar to previous work with these types of Coulomb clusters^, 
though three-dimensional systems have also been investigated using molecular dynamics 
simulations^. We also limit our investigation to a one-component plasma where particles 
have identical mass and charge. Other situations where the charges are not uniform have 
been investigated elsewhereP 

Our work can be stated in the following theorem: 

Theorem 1 For N charged particles, confined by a radially symmetric potential such that 
their equilibrium positions lie on a circle, the eigenfrequencies of the normal modes of oscil- 
lation for small perturbations about the equilibrium positions can be written as 

N 



r=2 



l,r) ( 1 — cos 



2n(h - l)(r - 1) 
N 



where §jv (hj) is 



, 3 + cos(f (i-j)) 



4^2 - 2 cos (f (i- j)) sin2 ~ 3)) ' 



In Section II we present the physical characteristics of the problem including our small- 
perturbation approximation in order to develop the conditions of Theorem]!} We develop the 
mass-weighted dynamical (Hessian) matrix based on the inter-particle Coulomb interaction 
using techniques which are typically used in normal mode analysis in physical chemistry^. 
We utilize concepts from graph theory to construct an edge-weighted Laplacian matrix 
whose weights correspond to the first-order spring constant for small oscillations about the 
equilibrium position. In Section III we apply the symmetry constraints on the equilibrium 
conditions and develop the properties of the edge weights Sjv (hj) for this physical situation. 
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We use these properties to prove that our dynamical matrix is a Laplacian matrix that is 
also symmetrically circulant. In Section IV we build on results from Gutierrez-Gutierre^^ 
to find the eigenvalues for a weighted Laplacian symmetrically circulant matrix. Finally, 
in Section V, we comment on the properties of the eigenvalues, the largest and smallest 
frequencies, and the limiting behavior as N becomes large. Since we limit the scope of this 
paper to the situation where the particle filling is low enough that there is only a single ring 
of particles with radius R in the equilibrium state, when we talk about larger numbers of 
particles, we assume that the trapping potential has been appropriately scaled to maintain 
this condition. 



II. PHYSICAL MODEL 



This paper will be focused on a system of identically charged particles confined to a two 
dimensional circle of radius R, as shown in Fig. [I] by an external potential. In order to 




FIG. 1: The equilibrium positions of ./V particles confined by an external potential to a circle of 
radius R. We identify each particle with an index i and define the angle 9i from the equilibrium 
position from the first particle. 

find the normal mode frequencies for small oscillations about this equilibrium position, we 
construct the dynamical (Hessian) matrix H based on the interaction Hamiltonian H of 
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this systemP The Hamiltonian has two pieces: the Coulombic repulsion between pairs of 
particles and a radially symmetric external potential. 

N N o N 



4vre \fi 



(1) 



8=1 



3=1 i>j 

It is natural to use polar coordinates, fl = (x^yi) = r^cos^, sin^j) because of the radially 
symmetric potential. The Hamiltonian in polar coordinates is thus 

2 N N ^ N 



H 



Q 

47re 



EE 

j=l i>j 



1 



rf + r 2 — 2riTj cos 



En 



(2) 



i=i 



We now assume that the radially symmetric external potential confines the particles to a 

Tj = R, and we obtain, 

2 -i - x ■ v / -i / , i , i , \ — 4 JV 



circle with the minimum at some radius R. Therefore 

N N 



H 



Q 



4vre 2R 



^EE 1 



j'=i *>i 



cos i 



i=l 



Using 



N 



V{R) = V(R) + V{R) + ... + V(R) = N (V(R)) 



(3) 



(4) 



i=l 

the interaction Hamiltonian is thus 

2 , JV N 



H 



Aixeq 2R 



EE 



1 — cos 



0; 



+ N(V(R)). 



(5) 



j=l i>j 

We can now construct the mass-weighted dynamical matrix consisting of the second partial 
derivatives of the Hamiltonian, evaluated at the equilibrium positions of the particles. This 
model is only valid for small perturbations about the equilibrium position and is only a 
first-order approximation. The mass weight makes the eigenvalues of the dynamical matrix 
the normal mode frequencies. Since we are working in polar coordinates, we rewrite the 
partial derivatives dxi Rd8i, and the elements of the dynamical matrix are 

1 d 2 H 



kl 



MR 2 de k d9i 

«2 



cquil 
d 2 



47re (2 J R 3 M) 86 h d9i 



' N N 

EE 



1 — cos (9j — 9a 



+ NV{R) 



(6) 



equil 



where k, I are integers and run from 1,...,N. We scale the dynamical matrix by the nat- 
ural frequency for this system, u 2 = g 2 /47re (2i? 3 M) which is equivalent to scaling the 



Hamiltonian by the harmonic oscillator frequency of the confining potential. Therefore all 
eigenfrequencies are scaled in these units as well. We also note that the NV(R) term is 
only dependent on R and so does not contribute after the partial differentiation. The scaled 
matrix elements are thus 



H, 



d 2 



ki 



d9 k d9i 



' N N 

EE 

.3=1 i>j 



1 - cos (0j - Oj) 



N N 

EE 



a 2 



3=1 i>j 



i d6 k d6i 



1 - cos (9i-9j) 



We can now simplify the partial derivatives using the Kronecker delta notation: 

{Sij - So). 



- 5 1 cot (-4^ 



d9i V 2 
We then take the second derivative, giving 



2 /l-cos(6» i -6» i ) 



(7) 



(8) 



d / 1 



cot 



89k \ 2 1 1-cosiOi-Qj) 



(Sij - S H ) 



3 + cos (6i - Oj) 



4 v /2-2cos(^-^) sin 2 f^ii 



(hj - Ski)(Sij - S u ). 



(9) 



Substituting this result into (|7| we obtain, 



N N 

ff« = EE 

j=l i>j 



3 + cos (9i — 9j] 



4 V / 2-2cos(^-^) sin 2,^— 



(SijSkj — 5ij5ki — 5n5kj + SuSki) 



J equil 



(10) 



In order to simplify our notation we will define the function § such that, 



Hi J) 



3 + cos (6i - 6j 



4 v /2-2cos(^-^) sin 2 f^i 



In graph theory § is the weighted edge function of the inter-connected particle inter- 
action graph. It is important to note that is an even function which leads to the 
following lemma: 

Lemma 1.1 The edge weight function S has the property that S = S (J,i)- 

Proof. For the cosine components of § (i, j), cos (9i — 9j) is equivalent to cos (— (Oj — Oi)) 
and because cosine is an even function is equivalent to cos (Oj — 9{). By symmetry this holds 
for the sin 2 component, because it is also even. Therefore S = S ■ 
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This also means that Hki = Hik, which agrees with the property that Hessian matrices are 
symmetric. This result is a consequence of Newton's third law since the interactions between 
two particles are equal. 

We will now define the edge weight function, S^v evaluated at the equilibrium posi- 
tions of the iV particles on the circle. The particles are equally distributed around the circle 
such that the difference in the angular position between any two particles is a multiple of 
2n/N . Therefore an arbitrary particle identified by the index i has an equilibrium angular 
position of 0j = 2iri/N + 6 Q where 8 is the arbitrary starting point on the circle (given by 
the angle of the first particle). Thus the difference in angular position between particles i 
and j is (0* - 6j) = 2tt (i - j) /N. 

Now that the equilibrium geometry is established we can define, the edge weight for N 
particles 

3 + cos(%(t-j)) 1 

§ JV («, J) = ■ 2/w YY" I 11 ) 



4^2 - 2 cos (f (z - j)) sin2 (U* - j)) ' 
An illustration of the first five interconnected particle graphs, connected by edges with the 



weight of S 7v (hj), is shown in Fig. [2} Using S^v (hj) we can further simplify (10). Beginning 

(1,2)\ 
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FIG. 2: The first five interconnected particle graphs showing one set of the weighted edges, §n (1, j) 
where j = 1, N. 



with the k I case and evaluating the Kronecker deltas, we are left with only two possibilities 

-S N (k,l) tfk>l 

(12) 

-§ N (l,k) i£k<l 

We can then combine this statement using Lemma |1.1[ Therefore the statement above can 



A7 



be written as 



Hki = —§at (k, I) for all k =^ I. 



(13) 
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Now we will look at the k = I case, from (10), 

( N 



fc-1 



Hkk — < 



J2s N (k,j) = ^(k,j) iik >3 

3=1 3=1 
N N 

^^§Af (i, k) = §(i,k) if k < i. 



(14) 



i=k+l 



we 



Recognizing the overlap between the two summations, and again using Lemma 1.1 
combine them to give 

N 

H kk = Yl MM (15) 
where the (j ^ k) notation means that the term j — k is skipped when performing the sum. 



Now (13) and (15) describe each entry of the matrix. 



III. PROPERTIES OF B N {i,j) 

Before we can prove that the weighted Hessian matrix generated by functions of 
is a Laplacian symmetric circulant matrix, we will first prove several properties of $n (hj)- 
We already know from Lemma |l.l[ that §^ = §jv CM) an d that Sjy is an even 
function. The new properties are contained in the following lemmas and their corollaries, 
with a graphical representation of the properties displayed in these proofs shown in Fig. [3j 

Lemma 1.2 The edge weight function §n has the property Sn (i, j) = $n (i + k, j + k) . 

Proof. The angle of the trig functions corresponding to Sat is 2tt (i — j) /N. Adding 
and subtracting a constant k from the argument gives us 2tt ((« + k) — (j + k)) /N, which is 
the angle corresponding to §jv (i + k,j + k). Therefore §n — §jv (i + k, j + k). g 

This result is a consequence that the interaction between any two particles separated by 
the same distance is the same (as shown in Fig. [3]) and corresponds to the discrete rotational 
symmetry of the physical system. 

Corollary 1.2.1 The edge weight function Sj\r(*>j) has the property $N{i,i + k) = 
§> N (i, i — k). 
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N-k+2 



FIG. 3: This diagram explicitly demonstrates Lemmas 1.2 and 3.2 The setup from which the 
other Lemmas can be graphically interpreted is also shown. 



Proof. By Lemma [l.l| §jy + k) is equivalent to §n {i + k,i). Using Lemma 1.2 



we can 



add an i — k to both sides, which yields Sjv (i, i — k). Therefore Sat (i, i + k) — §n — k). 



Lemma 1.3 The edge weight function §n (i,j) has the property §>n (0, j) = §n 

Proof. The angle of the trig functions corresponding to $n (0, j) is — 2tc (j) /N. Also the 
angle corresponding to Sn is — 2n (N — j) /N = 2n — (27ij/N). Therefore because the 

trigonometric components of Sjv (hi) are 2n periodic, S^v (0, j) = Sn (N,j). | 

There is no th particle on the circle; the physical interpretation of the Lemma is that one 
complete journey around the circle is achieved by a rotation of iV particles. For instance, 
starting at particle 1 on Fig. [3] and moving around N particles forward, brings one back to 
particle 1. 

Corollary 1.3.1 The edge weight function §Ar(i,j) has the property S^(i,k) = 
§ N (i,k±N). 



Proof. The function § N (i, k — N) by Corollary 1.2.1 is equivalent to S^r (i, k + N). Adding 



—k to both sides, by Lemma 1.2 we obtain Sjy (i — k,N). Which, by Lemma 1.3, is equal 



to Sat (i — k, 0). Adding k back to both sides yields Sat (i, k). 
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IV. LAPLACIAN SYMMETRIC CIRCULANT MATRIX 



Our analytic solution requires that the weighted Hessian H be a symmetric circulant 
matrix. We gain other insights about H if it is Laplacian. In this section we prove that the 
matrix is a Laplacian symmetric circulant matrix. 

Theorem 2 The matrix H is a Laplacian matrix, meaning that the sum of each row and 
column is zero. 

Proof. Because H is symmetric it suffices to show that the sum of each row is zero. The 
sum of each row can be written as 

N 

1=1 

Taking the I = k term out of the sum we obtain 

N 

Hkk + /J Hkh 



We can now substitute in the results from Eqns. (13) and (15); the sum of the row becomes 



N 



N 



N 



N 



MM) + Yl (-MM))= MM)- Y MM) = o. 

Therefore the sum of each row is zero. Thus the matrix H is a Laplacian matrix. | 

The properties of the eigenvalues of a Laplacian matrix are such that all of the eigenvalues 
are real and positive and at least one eigenvalue will be zero. The zero eigenfrequency for 
the normal mode oscillations corresponds to the common rotational mode of motion of the 
particles moving in the same direction around the circle.^ 



Theorem 3 The matrix H is symmetrically circulant, given by the for, 



circ7v(c ,ci, 



. Cjv-l) 



C C N _i Ctv_ 2 ••• Ci 

Ci Co C/v-i ' ■ : 

C2 Ci Co ' ■ Ctv-2 

: '•• '•• '•• cjv-i 

Cat_i ... C 2 Ci c 
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A symmetric circulant matrix matrix, B n , has the form that 



B„ 



circAr(&o, bi, ■ • ■ , bn-i , bn-i , ■ ■ ■ , 61) if n is odd 
circ N (b , 61, • • ■ ,6« i,6»,6» ,61) i/n zs even. 



To prove that matrix, H is a symmetric circulant matrix we need only prove that all of 
its diagonals are the same and the rows and columns follow the pattern for both odd and 
even cases. Physically we know that because all of the particles are identical the diagonals 
of the matrix should be the same. 

Lemma 3.1 The matrix H is a Toeplitz matrix, therefore the entries of each left-to-right 
diagonal is uniform. 

Proof. This proof will cover the different two cases for H k f. k 7^ I and k = I. 

For k 7^ I we know that from Lemma 1.2 S^v {hj) = §tv (i + k, j + k). Therefore, H k i 



is equal to H, t j where % = k + m and j = I + m for some integer m. Since a left-to-right 
diagonal is described by stepping both indices by the same number, the diagonal entries in 
H are the same for k 7^ I. 

For k — I, we will show that an arbitrary entry, H^k on the diagonal is equal to the next 
entry down the diagonal, H^ k+1 ^ k+1 y From Eq. (15) we know that the central diagonal is of 
the form 

N 

H kk = ^N(k,j) = S N (k,l) + S N (k,2) + ... + S N (k,N-l) + § N (k,N). (16) 

where k 7^ j. Also, the next entry down the diagonal is similarly 

N 

#(fe+i)(fc+i) = &N(k + l,j) (17) 

= §jv (k + 1, 1) + § N (k + 1, 2) + ... + Sat (k + 1, N - 1) + § N {k + 1, N) . 



Using Lemma L2 we can add —1 to j and k, giving 

N 

#(fc+i)(fc+i) = ^2 &N(k,j-l) 

= % N (k, 0) + S N (k, 1) + ... + S N (k, N - 2) + S N (k, N-l). (18) 
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All of the terms between Hub and H, 



f (fc+i)(fc+i) 



from (16) and (18) overlap except §N{k,N) 



from (16) and §n (k,0) from (18). However from Lemma 1.3, Sat (k,0) = Sjy (k,N). There- 
fore H k k = #(fc+i)(fc+i), and all of the central diagonals are the same. g 

In order to show that the matrix is symmetrically circulant, we also need the following 
lemma: 

Lemma 3.2 The first row of the matrix H has the symmetry such that $^(l,k) = 
S N (1,N- k + 2). 



Proof. By Corollary 1.3.1 Sat (1, iV — k + 2) is equivalent to Sat (1, —k + 2). Using Lemma 
1.2| we can add k — 1 to both sides, which yields §n (k, 1), which by Lemma [l.l| is equal to 
Stv (1, k). Therefore S^v (1, k) = S N (1, N - k + 2). B 

This property is demonstrated as the connection between the first particle and the k th 
particle going forwards and backwards around the circle in Fig. |3j 

With this proof, we have shown that the dynamical matrix H is symmetrically circulant 
of the form given in Theorem [3] g 

We now use the results by Gutierrez-Gutierrez and extend those results to find the eigen- 
values for a Laplacian symmetrically circulant matrix. The eigenvalues io\ (h = 1, . . . , N) 
of a symmetrically circulant matrix arep^ 



N+l 

+ H n cos ( 

r=2 ^ 

N 

H u + 2^2 H n cos ( 

r=2 ^ 



2n(h - l)(r - 1) 
N 

2n{h- l)(r- 1) 
N 



if is odd 

cos(7r(/i — 1)) if iV is even. 

(19) 



These two cases can be simplified to a single case by recognizing that the sum is doubled. 
Since the matrix has the form shown in Theorem [3j the sums in Eq. (19) are equivalent to 
a single sum over the entire row. There is no longer a distinction between the even and odd 
cases, giving a single sum: 

N 



rr V-rr /2tt(/i - 1) (r - 1) \ 
H u + Hrl COS — * i 

r=2 ^ ' 



(20) 
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We use our results for the explicit form of the dynamical matrix, Eqs. (13) and (15), and 
that H is symmetric, and find that 

^ = £ S * (1, r) - £ Siv (1, r) cos — ^ ^ ^ . (21) 

We then combine the two sums, since they sum over the same index and range. We also 
note that the eigenfrequencies have the same structure as the first row of the matrix. There 
are thus only \{N + l)/2] unique eigenfrequencies, where the notation |~] denotes the ceil- 
ing function, meaning that the number is rounded up to the nearest integer. The unique 
eigenfrequencies are thus 

N 



r=2 



E 



>N 



1, r) [ 1 — cos 



2?r(/i- l)(r - 1) 
N 



for h < \{N + l)/2]. (22) 



V. EIGENFREQUENCIES 



Finally, we evaluate the normal mode eigenfrequencies from Eq. ( 22 ) for large numbers of 
particles. Since there are approximately N/2 eigenvalues for iV particles, we evaluated the 
eigenfrequencies on a logarithmic scale, shown for the first 10,000 particles in Fig.|4} We have 
omitted the zero eigenfrequencies, as noted above, since that corresponds to the common 
rotation of all the particles around the circle. There is a clear frequency gap between the 
second and third frequency which means that, experimentally, these will always be uniquely 
addressable. However, spacing between the higher mode frequencies becomes smaller as the 
number of particles increases. 

The highest frequency mode is also of interest and we show a linear fit to this mode 
in the figure. The limiting behavior of this frequency is dominated by the denominator of 



the expression in (22). This denominator can be rewritten as the absolute value of a sine 
cubed, giving rise to an iV 3 dependence of the frequency. The solid line is a fit to u 2 = £N 3 
with a fit parameter of £ = 0.067838. Although the analytic evaluation of the highest mode 



frequency (h — > N/2 in Eq. (22)) is not trivial, is is possible to show that the fit parameter 
approaches 

g>> (23) 
as iV — > oo where ((3) is the Riemann Zeta function. 
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FIG. 4: Eigenfrequencies for the first 10,000 particles, as calculated from Eq. (22). The solid line 



is a fit to iV 3 for the highest frequency mode for large numbers of particles (N S> 1) 
VI. CONCLUSION 



We have proven Theorem 1, for N number of charged particles confined by a cylindrically 
symmetric potential to a circle. Small perturbations about the equilibrium positions give 
rise to normal modes of oscillation which are, 



where S/v is 



N 

E 

r=2 



'TV 



l,r) ( 1 — cos 



2tt(/i - l)(r 
N 



>7V {1,3, 



3 + cos(f (i-j)) 



4^2 - 2 cos (f - j)) Sin2 ~ j) ) ' 
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